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Abstract. This paper summarizes the research on the optimization 
and guidance of flight trajectories in the presence of windshear, 
performed by the Aero-Astronautics Group of Rice University 
during the period 1984-89. This research refers to windshear 
recovery systems and covers two major areas of investigation: 
optimal trajectories for take-off, abort landing, and penetration 
landing; and guidance schemes for take-off, abort landing, and 
penetration landing. 

Key Words. Flight mechanics, windshear problems, take-off, 
abort landing, penetration landing, optimal trajectories, optimal 
control, feedback control, guidance strategies, piloting strategies. 
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1 . Introduction 

Low-altitude windshear is a threat to the safety of 
aircraft in take-off and landing. Over the past 20 years, some 
30 aircraft accidents have been attributed to windshear. The 
most notorious ones are the crash of Eastern Airlines Flight 
066 at JFK International Airport (1975) , the crash of PANAM Flight 
759 at New Orleans International Airport (1982), and the 
crash of Delta Airlines Flight 191 at Dallas-Fort Worth 
International Airport (1985) . These crashes involved the loss 
of some 400 people and an insurance settlement in excess of 500 
million dollars. 

To offset the windshear threat, there are two basic systems: 
windshear avoidance systems and windshear recovery systems. A 
windshear avoidance system is designed to alert the pilot to 
the fact that a windshear encounter might take place; here, the 
intent is the avoidance of a microburst. A windshear recovery 
system is designed to guide the pilot in the course of a 
windshear encounter; here, the intent is to fly smartly across 
a microburst, if an inadvertent encounter takes place. Obviously, 
windshear avoidance systems and windshear recovery systems are 
not mutually exclusive, but complementary to one other. 

Examples of windshear avoidance systems are: ground-based 
mechanical systems (anemometers) , ground-based radar systems 
(Doppler radar), and airborne systems (radar or lidar). Examples 
of windshear recovery systems are: maximum angle of attack 
guidance, constant pitch guidance, acceleration guidance, and 
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gamma guidance. At this time, some of the above avoidance 
systems and recovery systems appear to be promising. Further 
research is both necessary and desirable in order to let the 
dust settle and prior to making large commitments of funds to one 
system or another. For previous research, see Refs. 1—117. 

1.1. Rice University Research on Windshear . This report 
refers to windshear recovery systems and summarizes the research 
performed at Rice University during the period 1984-89 under 
the sponsorship of NASA Langley Research Center, Boeing Commercial 
Airplane Company, and Air Line Pilots Association. This research 
was initiated in 1983 at the suggestion of Captain W. W. Melvin 
of Delta Airlines and ALPA. Its objective was: to study 
three problem areas, namely, take-off, abort landing, and 
penetration landing; for each problem area, to determine 
optimal trajectories, namely, trajectories minimizing a suitable 
performance index; for each problem area, to develop guidance 
schemes approximating the optimal trajectories in real time. 

While the above problem areas are quite different from 
one another, they are related by common mathematical grounds 
and common procedures. The optimization problems can be solved 
by means of a single algorithm, the sequential gradient-restoration 
algorithm, developed by the Aero-Astronautics Group of Rice 
University over the years 1970—85. From a comprehensive study 
of the optimal trajectories, the dominant properties of these 
trajectories can be found. Then, these dominant properties are 
employed to develop guidance laws that it is desirable to 
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approach in actual flight. Finally, these guidance laws are 
implemented via feedback control schemes in such a way that 
the guidance trajectories approximate the optimal trajectories. 

1.2. Wind Model . Although no two windshear encounters 
are exactly alike, two basic phenomena are always present: 
shear and downdraft. Therefore, it is important that these 
essential characteristics be present in the wind model employed 
in optimization and guidance studies. In this report, the 
assumed wind model has the following properties: (a) there is 

a transition from a uniform headwind to a uniform tailwind, 
with nearly constant shear in the core of the downburst; (b) 
the downdraft achieves maximum negative value at the center of 
the downburst; (c) the downdraft vanishes on the ground, h = 0; 
and (d) the wind velocity components nearly satisfy the continuity 
equation and the irrotationality condition in the core of the 
downburst . 

In this model, the horizontal shear (hence, the horizontal 
wind component) is independent of the altitude; and the down- 
draft (the vertical wind component) increases linearly with the 
altitude. Therefore, the wind model has the form 

W x = XA(x), W h = A(h/h*)B(x). (1) 

Here, the parameter A = AW /AW characterizes the intensity 

x x* 

of the windshear/downdraf t combination; the function A(x) 
represents the distribution of the horizontal wind versus the 
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horizontal distance; the function B(x) represents the distribution 
of the vertical wind versus the horizontal distance; and h* 
is a reference altitude, h* = 1000 ft. Also, AW x is the 
horizontal wind velocity difference (maximum tailwind minus 
maximum headwind) and AW = 100 fps is a reference value for 

‘ At 

AW . 
x 

Decreasing values of X (hence, decreasing values of AW x ) 
correspond to milder windshears; conversely, increasing values 
of X (hence, increasing values of AW^) correspond to more 
severe windshears. Therefore, by changing the value of X, 
one can generate shear/downdraft combinations ranging from 
extremely mild to extremely severe. 

Figure 1 shows the functions W^(x) and W^(x,h), computed 
for X = 1 , corresponding to AW = 100 fps. For other values of 
X (hence, for other values of AW^) , the ordinates in Fig. 1 
must be scaled proportionally to X (hence, proportionally to 
AW x ). 

1.3. Outline. The take-off problem is discussed in 
Section 2; the abort landing problem is considered in Section 3; 
and the penetration landing problem is presented in Section 4. 
Future developments are discussed in Section 5. Two appendices 
contain the equations of motion (Section 6) and the nomenclature 
(Section 7) . The paper ends with an extensive bibliography 
on the windshear problem; see Refs. 1-65 for research done 
at Rice University under the present grant, and see Refs. 66-117 
for research done elsewhere. 
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2. Take-Off Problem (Refs. 49-52, 54, 56, 61-65) 

2.1. Optimal Trajectories . Optimal take-off trajectories 
were studied with the aid of the sequential gradient-restoration 
algorithm for optimal control problems. Several performance 
indexes were considered; the most reliable one was found to be 
the deviation of the absolute path inclination from the nominal 
value. Because the deviation has a maximum value along the 
flight path, attention was focused on minimizing the peak 
deviation . The resulting optimization problem is a minimax 
problem or Chebyshev problem of optimal control, in which the 
desired minimum performance index has the form 

I* = min I, (2a) 

a 


I = max | Y 0 - Y e0 I , 

y e = arctan [ (VsinY + W h )/(VcosY + W x ) ] . 


(2b) 

(2c) 


Here, t is the running time; t is the final time; a is the 
angle of attack; Y e is the absolute path inclination; Y e0 
is the value at Y e at t = 0; Y is the relative path inclination; 
V is the relative velocity; and W , are the wind components 
[see Eqs . (1) ] . 

This problem was studied under the assumption that maximum 
power setting is employed and that inequality constraints are 
imposed on the angle of attack and the angle of attack rate. 
Hence, the power setting 3 is such that 


3 = 1 , 


3=0, 


(3a) 
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and the angle of attack a is such that 

« • * 

a <_ a*, -a* £ a _< a*, (3b) 

where a* is the stick-shaker angle of attack and a* = 3.0 deg/sec 
is the limiting angle of attack rate. 

With the above understanding, optimal trajectories 
were computed for three Boeing aircraft (B-727, B-737, and 
B-747) and several windshear intensities. From the extensive 
computations, certain general conclusions became apparent: 

(i) the optimal trajectories achieve minimum velocity 
at the end of the shear; 

(ii) the optimal trajectories require an initial decrease 
in the angle of attack, followed by a gradual increase; the 
maximum permissible angle of attack (stick-shaker angle of attack) 
is achieved near the end of the shear; 

(iii) for weak-to-moderate windshears, the optimal trajectories 
are characterized by a continuous climb; the average value of 

the path inclination decreases as the intensity of the shear 
increases ; 

(iv) for relatively severe windshears, the optimal 
trajectories are characterized by an initial climb, followed by 
nearly -horizontal flight, followed by renewed climbing after 
the aircraft has passed through the shear region; 

(v) weak-to-moderate windshears and relatively severe 
windshears are survivable employing an optimized flight strategy; 
however, extremely severe windshears are not survivable, even 
employing an optimized flight strategy; 
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(vi) in relatively severe windshears, optimal trajectories 
have a better survival capability than maximum angle of 
attack trajectories and constant pitch trajectories. 

For a particular case, the Boeing B-727 aircraft with a 
take-off weight W = 180000 lb, a flap deflection = 15 deg, 

and an initial altitude h^ = 50 ft, optimal trajectories were 
computed for 

A = 0.8, 1.0, 1.1, (4a) 

corresponding to 

AW = 80, 100, 110 fps. (4b) 

Figure 2 shows the resulting altitude profile h(t), velocity 
profile V ( t ) , and angle of attack profile a(t). Clearly, as 
AW^ increases, the altitude profile of the optimal trajectory 
changes: it is entirely ascending for AW^ = 80 fps; it includes 

a nearly-horizontal branch for AW^ = 100 fps; and it includes a 
descending branch for AW^ = 110 fps. Should the horizontal wind 
velocity difference be further increased to AW = 120 fps, 
the B-727 would crash, even flying an optimal take-off trajectory. 

2.2. Guidance Trajectories . The computation of the optimal 
trajectories requires global information on the wind field; 
that is, it requires the knowledge of the wind components at 
every point of the region of space in which the aircraft is 
flying. In practice, global information is not available; even 
if it were available, there would not be enough computing 
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capability onboard and enough time to process it adequately. 

As a consequence, the optimal trajectories are merely benchmark 
trajectories that it is desirable to approach in actual flight. 

Since global information is not available, what 
one can do is to employ local information on the wind 
field, in particular, local information on the wind acceleration 
and the downdraft. Therefore, the guidance problem must be 
addressed in these terms: Assuming that local information is 
available on the wind acceleration, the downdraft, as well as 
the state of the aircraft, we wish to guide an aircraft 
automatically or semiautomatically in such a way that the key 
properties of the optimal trajectories are preserved. 

Based on the idea of preserving the properties of the optimal 
tra jectories , three guidance schemes were developed at Rice 
University: (a) acceleration guidance, based on the relative 

acceleration; (b) absolute gamma guidance, based on the absolute 
path inclination; and (c) theta guidance, based on the 
pitch attitude angle. Among these, the acceleration guidance 
and the gamma guidance are of particular interest. 

Key to all of the guidance schemes is the recognition of 
the fact that, in the shear region, the dynamical effects due 
to shear and downdraft can be combined into a single factor, 
called the shear/downdraft factor (Ref. 54), 

F = w /g - w, /v. 

x h 


( 5 ) 
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The larger the value of F, the larger are the dynamical effects 
caused by the shear/downdraft combination. 

In the acceleration guidance (AG) , the guidance law 
has the form 

(AG) V/g + C F =0; (6) 

therefore, it enforces the proportionality between the shear/ 
downdraft factor and the instantaneous acceleration. In the 
gamma guidance (GG) , the guidance law has the form 

(GG) y e - Y e0 d - C 2 F) = 0; (7) 

therefore, it enforces the achievement of decreasing values 
of the absolute path inclination with increasing values of the 
shear/downdraft factor. Note that (6) and (7) apply to the 
shear region. 

In the aftershear region, guidance laws different from 
(6) and (7) are needed. Specifically, in the acceleration guidance 
(AG) , partial velocity recovery is enforced, 

(AG) V - C-V = 0. (8) 

In the gamma guidance (GG) , partial path inclination recovery 
is enforced, 

(GG) Y e - C 4 y e0 = 0. (9) 

These guidance laws are then implemented via feedback 
control forms, such that the difference a - a is made proportional. 
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via suitable gain coefficients, to the violation of any of 
the guidance laws ( 6 ) - ( 9 ) . Here, a is the instantaneous 
angle of attack and a = a (V) is the nominal angle of attack. 

Guidance trajectories were computed for the Boeing B-727 

aircraft using the above guidance schemes. It was found that 

both the acceleration guidance (AG) and the gamma guidance (GG) 

produce trajectories which are close to the optimal trajectories 

(OT) . In addition, the resulting near-optimal trajectories 

are superior to the trajectories arising from maximum angle of 

attack guidance (MAAG) and constant pitch guidance (CPG) . See 

Fig. 3, which refers to the Boeing B-727 aircraft,! = 1.0, 

and AW = 100 f ps . While the MAAG trajectory crashes, the CPG 
x 

trajectory scrapes the ground; on the other hand, the AG 
trajectory is close to the OT . An analogous remark holds for 
the GG trajectory, which is not shown, since it is nearly identical 
with the AG trajectory. 

2.3. Simplified Guidance Trajectories . As stated above, 
the previous guidance schemes require local information on the 
wind acceleration, the downdraft, and the state of the aircraft. 
While this information will be available in future aircraft, it 
might not be available on current aircraft. 

For current aircraft, one way to survive a windshear encounter 
is the quick horizontal flight transition technique, based on 
the properties of the optimal trajectories. The quick horizontal 
flight transition technique requires an initial decrease of 
the angle of attack, so as to decrease the path inclination to 
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nearly horizontal. Then, near ly-horizontal flight is maintained 
during the windshear encounter. Climbing flight is resumed 
after the shear is past. 

For relatively severe windshears, the quick horizontal 

flight transition technique yields trajectories which are 

competitive with those of the guidance schemes discussed previously. 

In addition, for relatively severe windshears, the quick 

horizontal flight transition technique yields trajectories 

which have better survival capability than those resulting 

from maximum angle of attack guidance and constant pitch guidance. 

2.4. Survival Capability . Perhaps, the best way of 

assessing the merit of a particular guidance scheme is the 

concept of survival capability. Consider the one parameter 

family of wind models (1) . As the parameter X increases, more 

intense windshear/downdraf t combinations are generated until a 

critical value of X is found (hence, a critical value of AW x 

is found), such that h . =0 for a given trajectory type. Thus, 

mm 

the survival capability is the critical wind velocity difference 

AW for which a crash first occurs, 
xc 

The results are shown in Table 1, which supplies the 

survival capability AW for the optimal trajectory and various 

X c 

guidance trajectories. Table 1 also shows the windshear efficiency 
ratio h , defined to be 


= (AW ) nm /(AW 

xc PT xc 


OT ’ 


n 


( 10 ) 
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Here, the subscript PT denotes a particular trajectory and the 
subscript OT denotes the optimal trajectory. Clearly, if the 
windshear efficiency of the OT is defined to be 100%, that of 
the AG trajectory is 95%, that of the GG trajectory is 96%, that 
of the CPG trajectory is 85%, and that of the MAAG trajectory 
is 48 % . 
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3. Abort Landing Problem (Refs. 53, 57-59, 61-65) 

3.1. Optimal Trajectories . Optimal abort landing 
trajectories were studied with the aid of the sequential 
gradient-restoration algorithm for optimal control problems. 

Several performance indexes were considered; the most reliable 
one was found to be the drop of altitude with respect to a 
reference value. Because the altitude drop has a maximum value along 
the flight path, attention was focused on minimizing the peak 
altitude drop. The resulting optimization problem is a minimax 
problem or Chebyshev problem of optimal control, in which the 
desired minimum performance index has the form 


I* = min I, 

a 

I = max | h R - h | , 


(11a) 

0 < t < t. (lib) 


Here, t is the running time; x is the final time; a is the 
angle of attack; h is the altitude above ground; and h R is a 
reference value for the altitude, for instance h R = h* = 1000 ft. 

This problem was studied under the assumption that, at the 
windshear onset, the power setting 6 is increased from the 
initial value 3 = Bq to the maximum value 3 = 1 at a constant 
time rate 3 = 3^; afterward, the maximum value is kept. Inequality 
constraints are imposed on the angle of attack and the angle 
of attack rate. Hence, 


e = e 0 + So 1 ' 


0 < t < a, 


(12a) 
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g = l CJ £ t £ x , (12b) 

with a = (1 - B 0 )/i3 0 , and 

a £ a*, -a* £ a 1 a*' (12c) 

where a* is the stick-shaker angle of attack and a* = 3.0 deg/sec 
is the limiting angle of attack rate. 

With the above understanding, optimal trajectories 
were computed for three Boeing aircraft (B-727, 

B-737, and B-747) and for several combinations of windshear 
intensity, initial altitude, and power setting rate. From the 
extensive computations, certain general conclusions became 
apparent with reference to strong-to-severe windshears : 

(i) the optimal trajectory includes three branches: a 
descending flight branch, followed by a nearly-horizontal 
flight branch, followed by an ascending flight branch after 
the aircraft has passed through the shear region; 

(ii) along an optimal trajectory, the point of minimum 

velocity is reached at the end of the shear; 

( iii) the peak altitude drop depends on the windshear 
intensity, the initial altitude, and the power setting rate; it 
increases as the windshear intensity increases and the initial 
altitude increases; and it decreases as the power setting 

rate increases; 

(iv) the peak altitude drop of the optimal 
trajectory is less than the peak altitude drop of the maximum 
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angle of attack trajectory and the constant pitch trajectory; 

(v) the survival capability of the optimal trajectory 
is superior to that of the maximum angle of attack trajectory 
and the constant pitch trajectory. 

For a particular case, the Boeing B-727 aircraft with a 
landing weight W = 150000 lb, a flap deflection = 30 deg, 

and an initial altitude h^ = 600 ft, optimal trajectories were 
computed for 


A = 1.0, 1.2, 1.4, 
corresponding to 


(13a) 


AW X - 100, 120, 140 f ps . (13b 

Figure 4 shows the resulting altitude profile h (t) , velocity 

profile V(t), and angle of attack profile a(t). Clearly, as 

AW x increases, the altitude profile of the optimal trajectory 

changes; in particular, as AW^ increases, the minimum altitude 

of the optimal trajectory decreases, even though the aircraft 

still survives the windshear encounter. Should the horizontal 

wind velocity difference be further increased to AW = 190 fns 

x ^ r 

the B-727 would crash, even flying an optimal abort landing 
trajectory. 

3 - 2 * Guidance Trajectories . By necessity, guidance 
trajectories must rely on only local information on the wind 
acceleration, the downdraft, and the state of the aircraft. 

The intent is to guide an aircraft automatically or semi- 
automa tic ally in such a way that the key properties of the 



16 


AAR- 2 4 4 


optimal trajectories are preserved. 

Based on the idea of preserving the properties of the 
optimal trajectories, five guidance schemes were developed at 
Rice University: (a) target altitude guidance, based on the 

initial altitude and the total wind velocity difference; (b) 
safe target altitude guidance, based on only the initial 
altitude; (c) acceleration guidance, based on the relative 
acceleration; (d) gamma guidance, based on the absolute path 
inclination; and (e) theta guidance, based 

on two target pitch angles, a lower target pitch followed by 
a higher target pitch. Among these, the acceleration guidance 
and the gamma guidance are of particular interest. 

Key to all of the guidance schemes is the recognition of 
the fact that, in the shear region, the dynamical effects due 
to shear and downdraft can be expressed via the shear/downdraft 

factor (Ref. 54) 

F = W /g - W, /V. (14) 

x h 

The larger the value of F, the larger are the dynamical effects 
caused by the shear/downdraft combination. 

In the acceleration guidance (AG) , the guidance law has 
the form (6) for the shear region and the form (8) for the 
aftershear region. For the shear region, the constant C 1 takes 
two values, one for the descending flight branch and one for 
the nearly-horizontal flight branch. Appropriate switch conditions 
regulate the transition from one branch to another of the AG 
trajectory . 
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In the gamma guidance (GG) , the guidance law has the form 
(7) for the shear region and the form (9) for the aftershear 
region. For the shear region, the constant C 2 takes two 
values, one for the descending flight branch and one for the 
nearly-horizontal flight branch. Appropriate switch conditions 
regulate the transition from one branch to another of the GG 
trajectory . 

These guidance laws are then implemented via feedback 
control forms, such that the difference a - a is made proportional, 
via suitable gain coefficients, to the violation of any of the 
guidance laws (6) -(9). Here, a is the instantaneous angle of 
attack and a = a (V) is the nominal angle of attack. 

Guidanc-e trajectories were computed for the Boeing B-727 

aircraft using the above guidance schemes. It was found that 

both the acceleration guidance (AG) and the gamma guidance (GG) 

produce trajectories which are close to the optimal trajectories 

(OT) . In addition, the resulting near-optimal trajectories are 

superior to the trajectories arising from maximum angle of 

attack guidance (MAAG) and constant pitch guidance (CPG) . See 

Fig. 5, which refers to the Boeing B-727 aircraft, X = 1.2, 

and AW = 120 f ps . While the MAAG trajectory crashes, the CPG 
x 

trajectory survives, albeit with a minimum altitude about half 
that of the OT. On the other hand, the AG trajectory is 
close to the OT. An analogous remark holds for the GG trajectory, 
which is not shown, since it is nearly identical with the AG 


trajectory . 
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3.3. Survival Capability . Perhaps, the best way of 
assessing the merit of a particular guidance scheme is the concept 
of survival capability. Consider the one-parameter family 
of wind models (1) . As the parameter X increases, more 
intense windshear/downdraf t combinations are generated until 
a critical value of X is found (hence, a critical value of 
AW x is found) , such that h^^ n = 0 for a given trajectory type. 
Thus, the survival capability is the critical wind velocity 
difference AW for which a crash first occurs. 

X.C 

The results are shown in Table 2, which supplies the 
survival capability AW for the optimal trajectory and various 
guidance trajectories. Table 2 also shows the windshear 
efficiency ratio n , defined to be 


n 


“ (AW xc>Pt/ ( ™xo>OT- 


(15) 


Here, the subscript PT denotes a particular trajectory and the 
subscript OT denotes the optimal trajectory. Clearly, if the 
windshear efficiency of the OT is defined to be 100%, that of 
the AG trajectory is 96%, that of the GG trajectory is 98%, 
that of the CPG trajectory is 75%, and that of the MAAG trajectory 
is 44%. 
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4. Penetration Landing Problem (Refs. 55, 60, 62, 65) 

4.1. Optimal Trajectories . Optimal penetration landing 
trajectories were studied with the aid of the sequential 
gradient-restoration algorithm for optimal control problems. 

The performance index being minimized is the least-square 
deviation of the flight trajectory from the nominal trajectory. 
The resulting optimization problem is a Bolza problem of optimal 
control, in which the desired minimum performance index has the 
form 


I 


* 


I 


min I , 

a, 3 

T ~ 2 
[h - h (x) ] dt . 

0 


(16a) 

(16b) 


Here, t is the running time; x is the final time; a is the 

angle of attack; 3 is the power setting; x is the horizontal 

distance; h is the altitude above ground; and h(x) is the 

nominal altitude. In turn, the function h(x) is computed by 

assuming that the nominal trajectory includes two parts: the 

approach part (h > 50 ft) , in which the slope is constant, 

y = -3.0 deq; and the flare part (h < 50 ft), in which the slope 
e — 

y is a linear function of the horizontal distance, 
e 

This problem was studied under the assumption that the 
absolute path inclination at touchdown is to be Y e = -0.5 deg; 
that the touchdown velocity is to be within + 30 knots from 
the nominal value; and that the touchdown distance is to be 
within + 1000 ft of the nominal value. Inequality constraints 
were imposed on the power setting, the power setting rate, the 
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angle of attack, and the angle of attack rate. Hence, 

B* £ B < 1, -B* < B < B*, (17a) 

and 

a<a„ (17b) 

Here, B* is the lower bound for the power setting; B* i s 

the limiting power setting rate; a* is the stick-shaker angle 

of attack; and a* = 3.0 deg/sec is the limiting angle of attack 

rate . 

With the above understanding, optimal trajectories were 
computed for the Boeing B-727 aircraft and for several combinations 
of windshear intensity, initial altitude, and power setting 
rate. From the extensive computations, certain general conclusions 
became apparent: 

(i) the angle of attack has an initial decrease, which 
is followed by a gradual, sustained increase; the largest value 
of the angle of attack is attained near the end of the shear; in 
the aftershear region, the angle of attack decreases gradually, 

(ii) initially, the power setting increases rapidly 
until maximum power setting is reached; then, maximum power 
setting is maintained in the shear region; in the aftershear 
region, the power setting decreases gradually; 

( iii) the relative velocity decreases in the shear region 
and increases in the aftershear region; the point of minimum 
velocity occurs at the end of the shear; 
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(iv) depending on the windshear intensity and the initial 
altitude, the deviations of the flight trajectory from the 
nominal trajectory can be considerable in the shear region; 
however, these deviations become small in the aftershear region, 
and the optimal flight trajectory recovers the nominal trajectory; 

(v) the optimal trajectory is better able to satisfy 
the touchdown requirements concerning the absolute path 
inclination, the velocity, and the distance than the fixed control 
trajectory (fixed angle of attack, coupled with fixed power 
setting) and the autoland trajectory (angle of attack controlled 
via path inclination signals, coupled with power setting 
controlled via velocity signals). 

For a particular case, the Boeing B-727 aircraft with a 
landing weight W = 150000 lb, a flap deflection 6^ = 30 deg, 
and an initial altitude h^ = 600 ft, optimal trajectories were 
computed for 

X = 1.0, 1.2, 1.4, (18a) 

corresponding to 

AW = 100, 120, 140 fps. (18b) 

X 

Figure 6 shows the resulting altitude profile h(x), velocity 
profile V(x), angle of attack profile a(x), and power setting 
profile 8 (x) . 

4.2. Quasi-Op timal Trajectories. Quasi-optimal trajectories 


can be generated by minimizing the functional (16) under the 
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assumption that only one control is available, the angle of 
attack a(t), which is subject to Ineqs. (17b). The power 
setting 3(t) is specified a priori, based on the results on 
the optimal trajectories. 

For the shear portion of the trajectory, maximum power 
setting is maintained; that is, the function B(t) is given by 

3 = 3 q + 3 Q t, 0 < t < a, (19a) 

3 = 1 , a £ t £ u), (19b) 

• 

with a = (1 - 3 q )/3q. Here, t = 0 is the initial time, 
t = o is the time at which maximum power setting is reached, 
and t = a) is the time at which the shear terminates; also, 

3q is the initial power setting and 3 q is the initial power 
setting rate. 

For the aftershear portion of the trajectory, maximum 
power setting is maintained if V < Vj, and reduced power setting 
is maintained if V > V., where V is the instantaneous relative 
velocity and V £ = V Q - 30 knots is the lower bound for the 
velocity. Hence, 3(t) is supplied by the following relationship 

3 - 3 0 = (1 - S Q ) (v 0 “ v >/(v 0 - V' (19c) 

3* £ 3 £ 1/ w <_ t £ x. (19d) 

For a particular case, the Boeing B-727 aircraft with a 
landing weight W = 150000 lb, a flap deflection 6^, = 30 deg, 
and an initial altitude h Q = 600 ft, quasi-optimal trajectories 
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were computed for the conditions (18) . From the numerical 
results, upon comparing the optimal trajectory (OT) , the 
quasi-optimal trajectory (QOT) , and the nominal trajectory (NT) , 
the following conclusions became apparent: 

(i) the OT and the QOT are geometrically close to one 

another ; 

(ii) both the OT and the QOT are geometrically close 

to the NT, providing the windshear is not exceptionally severe, 

AW < 120 fps; 
x — 

(iii) both the OT and the QOT satisfy the touchdown 
requirements concerning the absolute path inclination, the 
velocity, and the distance. 

Results (i)-(iii) imply that, for the purposes of constructing 
a guidance scheme, the coupling relation between the angle of 
attack and the power setting can be ignored. This separation 
result simplifies to a considerable degree the design of 
guidance and control systems capable of approximating the 
behavior of the optimal penetration landing trajectory in the 
presence of windshear (Ref. 60). 

4 . 3 . Guidance Trajectories . A penetration landing 
guidance (PLG) scheme was constructed, based on the following 
ideas: to rely on only local information on the wind acceleration, 
the downdraft, and the state of the aircraft; to preserve the 
key properties of the optimal trajectories; to treat the power 
setting and the angle of attack as decoupled controls. Because 
of the separation result of the previous section, the power 
setting determination can be based on the velocity, while the angle 
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of attack determination can be based on the absolute path 
inclination and the glide slope. 

For the power setting, the guidance law has the form 

6-1=0, (20a) 

in the shear portion of the trajectory. This means that the 
power setting must be increased as soon as possible to the 
maximum value. 

Also for the power setting, the guidance law has the form 

V - V Q = 0, (20b) 

in the aftershear portion of the trajectory. This forces the 
achievement of velocities consistent with the velocity touchdown 
requirement . 

The guidance laws (20) are then implemented via a single 
feedback control form, such that the difference 6 - is a 
linear combination, via suitable gain coefficients, of two 
terms. The first term is proportional to the shear/downdraf t 
factor F and takes care indirectly of (20a) . The second term 
is proportional to the violation of the guidance law (20b) ; 
indeed, it is identical with the right-hand side of (19c) . 

For the angle of attack, the objective of the guidance law 
is to maintain the flight trajectory close to the nominal 
trajectory via two signals: the absolute path inclination Y e 
and the glide slope y . 

For the approach part (h 50ft) , the glide slope signals 
are predominant. Hence, the guidance law is dominated by 
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v - v (h , F) = 0. (21a) 

g g 

Here, y is the instantaneous glide slope; y is the nominal 
g 9 

value of y ; h is the instantaneous altitude; and F is the 

g 

shear/downdraft factor. The function Y a (h,F) is determined 
from the study of the optimal trajectories; for details, see 
Ref. 60. Enforcement of (21a) is essential to achieving 
the required touchdown distance. 

For the flare part (h < 50 ft) , the absolute path inclination 
signals are predominant. Hence, the guidance law is dominated 
by 

Y - Y (h,F) = 0. ( 21b ) 

e e 

Here, Y e is the instantaneous absolute path inclination; y g 
is the nominal value of Y e ? h is the instantaneous altitude; 
and F is the shear/downdraft factor. The function Y e (h,F) 
is determined from the study of the optimal trajectories; for 
details, see Ref. 60. Enforcement of (21b) is essential 
to achieving the required touchdown absolute path inclination. 

The guidance laws (21) are then implemented via a feedback 
control form, such that the difference a - a is a linear 
combination, via suitable gain coefficients, of the violations 
of the guidance laws (21). Here, a is the instantaneous angle 
of attack and ct = 5(V) is the nominal angle of attack. 

Guidance trajectories were computed for the Boeing B-727 
aircraft using the above guidance scheme. It was found that 
the penetration landing guidance (PLG) produces trajectories 
which are relatively close to the optimal trajectories (OT) . 
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In addition, the resulting near-optimal trajectories are superior 
to the trajectories arising via fixed-control guidance (FCG) 
and autoland guidance (ALG) in terms of the ability to meet 
the touchdown requirements and in terms of survival capability 
in strong-to-severe windshears. See Fig. 7, which refers to 
the Boeing B-727 aircraft, X = 1.2, and AW^ = 120 fps. Both 
the FCG trajectory and the ALG trajectory crash; however, the 
PLG trajectory survives the windshear encounter, just as the OT. 

4.4. Simplified Guidance Trajectories . From a practical 
point of view, it must be emphasized that penetration landing 
makes sense only if the windshear encounter occurs at lower 
altitudes; if the windshear encounter occurs at higher altitudes, 
abort landing must be preferred. Therefore, low-altitude 
penetration landing deserves particular attention. 

For this special situation, h Q <_ 200 ft, the guidance 
scheme of the previous section can be simplified by keeping 
the power setting at the maximum permissible value and by 
controlling the angle of attack via absolute path inclination 
signals. This is the same as stating that, among all the touch- 
down requirements, priority must be given to the touchdown 
path inclination requirement. 

For the power setting, the guidance law has the form 

6 - 1 = 0 , ( 22 ) 

which is now employed in both the shear portion and the aftershear 
portion of the trajectory. This means that the power setting 
must be increased as soon as possible to the maximum value and 
then kept at the maximum value until touchdown. 



27 


AAR- 2 4 4 


For the angle of attack, the objective of the guidance 
law is to maintain the flight trajectory close to the nominal 
trajectory via only path inclination signals. The guidance 
law has the form 

Y e " Y e (h) = 0, (23) 

which differs from (21b) in that the nominal absolute path 
inclination does not involve the shear/downdraft factor. There- 
fore, the function y (h) in (23) is the same as the function 
characterizing the nominal trajectory (NT) ; for details, see 
Ref. 6 0. 

The guidance law (23) is then implemented via a feedback 
control form, such that the difference a — ot is made proportional, 
via a suitable gain coefficient, to the violation of the guidance 
law (23) . Here, a is the instantaneous angle of attack and 
a = a(V) is the nominal angle of attack. 

Simplified guidance trajectories were computed for the 
Boeing B-727 aircraft using the above guidance scheme. It was 
found that the simplified penetration landing guidance (SPLG) 
produces trajectories which are close to both the optimal 
trajectories (OT) and the nominal trajectories (NT) . For 
detailed results, see Ref. 60.. 
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5. View Toward the Future 

Over the past five years, considerable research has been 
performed at Rice University on two aspects of the windshear 
problem: determination of optimal trajectories and development 
of near-optimal guidance schemes. It now appears that, over 
the next 2-3 years, an advanced windshear control system can 
be developed, capable of functioning in different wind models 
and capable of covering the entire spectrum of flight conditions, 
including take-off, abort landing, and penetration landing. 

The advanced windshear control system must combine an 
advanced windshear detection system and an advanced windshear 
recovery system. The advanced windshear detection system requires 
the use of real-time identification techniques and must be 
characterized by small computational time, coupled with limited 
memory requirements. The advanced windshear recovery system 
must incorporate four basic properties explained below: completeness, 
continuation, near-optimality, and simplicity. 

Completeness means that the system should be able to function 
in take-off, abort landing, and penetration landing. Continuation 
means that the system should cover a variety of situations, 
ranging from zero windshear to moderate windshear to strong— to— 
severe windshear; the switch from no— windshear operation to 
windshear operation should be smooth. Near— optimality means 
that the system should be constructed so as to supply a good 
approximation to the properties of the optimal trajectories. 
Simplicity means that the system should be as simple as possible 
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and should emphasize the use of existing instrumentation, 
whenever possible. 

With reference to the continuation property, it must be 
noted that any adverse wind gradient is both preceded and 
followed by a favorable wind gradient. The advanced windshear 
recovery system must not only react in a near-optimal way to 
adverse wind gradients, but must exploit to the best advantage 
of the aircraft favorable wind gradients. This means that, 
in an increasing headwind scenario, kinetic energy must be 
increased; conversely, in a decreasing tailwind scenario, 
potential energy must be increased. Clearly, this requires 
that not only the current windshear signals be measured, but 
that previous windshear signals be recorded and memorized, 
such that favorable wind gradients can be detected and utilized. 
Hence, some modification of the guidance schemes described in 
Sections 2-4 is in order. 

To sum up, it is felt that an advanced windshear control 
system, endowed with the properties described above, should 
improve considerably the survival capability of the aircraft 
in a severe windshear. 
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6 . Appendix A; Equations of Motion 

In this report, we make use of the relative wind-axes 
system in connection with the following assumptions: (a) the 

aircraft is a particle of constant mass; (b) flight takes place 
in a vertical plane; (c) Newton's law is valid in an Earth-fixed 
system; and (d) the wind flow field is steady. 

With the above premises, the equations of motion include 
the kinematical equations 

x = Vcosy + W x , (24a) 

h = Vsiny + W (24b) 

and the dynamical equations 

V = (T/m)cos(a + 6) - D/m - gsiny 

» • 

- (W x cosy + W^siny ) , (25a) 

y = (T/mV)sin(a +6) + L/mV - (g/V)cosy 

• • 

+ ( 1/V) (W x siny - Wxcosy) . (25b) 

Because of assumption (d) , the total derivatives of the wind 
velocity components and the corresponding partial derivatives 
satisfy the relations 

• 

w = ow /3x) (Vcosy + W ) + (9W/8h) (Vsiny + W.) , (26a) 

XX XX il 

W h = (3W h/ 3x) ( Vcos Y + w > + ( 9W h /8h) (Vsiny + W h ) . (26b) 
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These equations must be supplemented by the functional 


relations 

T = T (h, V, 0) , (27a) 

D = D (h , V, a) , L = L (h,V,a) , (27b) 

W x =W x (x,h), W h =W h (x,h), (27c) 

and by the analytical relations 

0 = a + y, (28a) 

y = arctan [ ( VsinY + W, )/(Vcosy + W )]- (28b) 

6 il X 


The differential system (24) -(27) involves four state 
variables [x(t), h(t), V(t), y ( t ) ] and two control variables 
[a(t), 8 (t) ] . However, the number of control variables reduces 

to one (the angle of attack) , if the power setting is specified 
in advance. The quantities (28) can be computed a posteriori, 
once the values of the state and the control are known. 
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7 . Appendix B: Notations 

Throughout the report, the following notations are employed 
D = drag force, lb; 

F = shear/downdraft factor; 

-2 

g = acceleration of gravity, ft sec ; 

h = altitude, ft; 

L = lift force, lb; 

-1 2 

m = mass, lb ft sec ; 

„ ^ . 2 

5 = reference surface area, ft ; 
t = running time, sec; 

T = thrust force, lb; 

V = relative velocity, ft sec 1 ; 

W = mg = weight, lb; 

= h-component of wind velocity, ft sec ; 

-1 

W^. = x-component of wind velocity, ft sec ; 
x = horizontal distance, ft; 
a = angle of attack (wing) , rad; 

3 = engine power setting; 

y = relative path inclination, rad; 

Y e = absolute path inclination, rad; 

= glide slope angle, rad; 

6 = thrust inclination, rad; 

6 = flap deflection, rad; 

0 = pitch attitude angle (wing) , rad; 

A = wind intensity parameter; 
t = final time, sec. 
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Table 1. Take-off survival capability for the B-727 aircraft, 
W = 180000 lb, 6 = 15 deg, h n = 50 ft. 

r U 


Trajectory 

A 

c 

AW 

xc 

(fps) 

n 

OT 

1.195 

119.5 

1.000 

AG 

1.130 

113.0 

0.946 

GG 

1.153 

115.3 

0 .965 

CPG 

1.018 

101.8 

0.852 

MAAG 

0.577 

57.7 

0 .483 

Table 2 . 

Abort landing survival capability for the 

W = 150000 lb, 5 = 30 deg, h_ - 600 ft. 

F l) 

B-727 aircraft. 

Trajectory 

A 

c 

AW 

xc 

(fps) 

n 

OT 

1.871 

187.1 

1.000 

AG 

1.791 

179.1 

0 . 957 

GG 

1.842 

184.2 

0.985 

CPG 

1.394 

139 . 4 

0.745 

MAAG 

0.817 

81.7 

0.437 
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